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SUMMARY 

An evaluation has been made of two flow-field codes employing the tiine- 
and space-marching numerical technicpies. The space -marching technique 
considered is ^ul implicit finite-difference method due to Davis. The 
time-marching method employs the time-asymptotic, two-step, finite-difference 
method of MacCormack. These two methods have been used to analyze the 
flow field around a massively blown Jupiter entry probe under perfect*gas 
conditions. In order to obtain a direct point-by -point comparison, the 
computations by the two schemes have been done by using identical grids 2md 
tiirbulence v>^els. For the same degree of accuracy, the space-marching scheiie 
takes much less time as compared to the time-marching method and would appear 
to provide acc\irate results for the problems with nonequilibrium chemistry, 
free from the effect of local differences in time on the final solution which 
is inherent in time-marching methods. With the time-marching method, however, 
the solutions are obtainable for the realistic entry probe shapes with massive 
or uniform surface blowing rates , whereas with the space -marching technique , 
it is difficult to obtain converged solutions for such flo\ conditions. The 
choice of the numerical method is, therefore, problem dependent. Both the 
methods give equally good results for the cases where a comparison with the 
experimental data has been made. 


*NRC-Senior Research Associate. 
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INTRODUCTION 

1 2 3 4 

Recently an evaluation of the three different flow-field codes ' ' 

was presented. Ihese codes were primarily developed for the hypersonic, 

massively blown, and strongly radiating shock layers with application to 

the Jupiter entry conditions. Evaluation and comparison of the codes are 

important because of the dependence of the Galileo heat-shield design on 

these flowr-field solutions. However, direct comparisons have been hankered 

by the fact that the three codes use different solution procedures, different 

cc»nputational mesh sizes, different chemical and turbulence models, different 

convergence criteria, as well as a different number of computational nodes 

along and perpendicular to ^e body. Even for the only cenmnon feature of 

radiative transport code, there are subtle differences in the spectral 

details of some species. Thus, for an objective evaluation of the different 

numerical solution methods employed by these codes, it would be desiredsle to 

select a simple no-biowing perfect-gas flow-field case for which the tvirbulent 

models are well established. The purpose of this paper is to present the 

results of such a study. Ihe results have been obtained by employing two of 

2 

the three numerical codes mentioned earlier. Ihe first of these two is a 
space-marching method that solves the steady-state viscous-shock-layer type 
equations by the method of Davis. ^ The second one^ is a time-marching method 
that employs the time-asymptotic , two-step, finite-difference method of 
MacCormack^ for the solution of time -dependent viscous-shock- layer type 

4 

equations. The third numerical code uses an inverse analytical method. 

Since thi method of reference 4 also falls under the category of a space-marchj.ng 
method i similar to the category of method of ref. 2) , it has not been 
included in the present compariscxi. The work presented here includes the 
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results for a blown perfect-gas viscous shock layer by the space- amd 
time-marching numerical methods. The comparative study of these two methods 
shows their relative merits and disadvantages for a blown and unblown shoc>’ 
layer. Ihe present results may also be helpful in resolving the differences^ 
in predicted results by the different codes for the Jupiter entry conditions. 

ANALYSIS 

2 

As mentioned earlier, the space-marching (perfect-gas version of HYVIS ) 
and the time-marching codes (perfect-gas version of COLTS^) have been employed 
to predict the perfect-gas forebody flow field for Jupiter entry-type conditions. 
For this study, these two codes employ the same grid spacing along and perpen- 
dicular to the body, and the same number of computational mesh nodes in the 
streamwise and normal directions. The flow for both of these methods under- 
goes transition instantaneously to turbulent flow at the first body station 

downstream of the stagnation point. The two methods also employ the same 
7 8 

turbulence model ' with and without blowing. For the case which includes 

blowing, identical blowing-rate distribution as shown in figure 1 has been used. 

The governing equations employed in the present study can be obtained 

from the unsteady/steady Navier-Stokes equations by keeping terms up to second 

order in the inverse square root of the Reynolds number in both the viscous 

and inviscid regions. These equations, when modified for turbulent flow by 

using methods analogous to those presented by Dorrance,^ provide perfect-gas 

8 10 

governing equations for the space- and time -marching methods. 

In addition to the surface blowing rate distribution of figure 1, the 
boundary conditions employed along the body surface are no velocity slip and 
temperature jump. The wall temperature is taken as a constant specified 
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value. For the tine-marching method, the wall pressure is taken equal to 
the pressure at the adjacent grid point in the normal direction. TTie 
conditions imposed at the shock are calculated by using the Rankine-Hugoniot 
shock relations. Ihese relations in the body-oriented coordinates (for use 
wi'h the time-marching method) are provided in reference 11. The flow conditions 
along the supersonic outflow boundary are determined by extrapolation frcxn 
the inner mesh points for the time-marching scheme. No such flow conditions 
are required for the space-Ti.«rc U\g method. Further, for the present study. 



has been used to determine the viscosity, and the coefficient of thermal 
conductivity is computed by assuming a constant Premdtl number: 

k* » y 

(y - l)Pr 

In equations (1) and (2) , the quantities with an asterisk are dimensional and R* 
is the gas constant. 

In order to obtain a direct point-by-point comparison and to keep the 
seune truncation error in the finite-difference approximations employed in 
the space- amd time -inarching schemes, the computations by the two schemes 
have been done by using identical grids. The following logarithmic 
grid-spacing 

in 

n » 

(3) 


6 ♦ (1 - 7 ) 

6 - 0 - 7 > 


in 
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has been used In the two mathods in the direction norroal to the body to 

allow higher resolution near the body surface. lo equation (3)# $ is the grid- 

• * * * 

stretching factor, n(*n /r„) is the normal coordinate, and 6(«6 /r„) is 

N N 

the shock standoff distance. Finally, the following convergency criterion 
has been employed for both blown as well as unblown cases: 


maximum 

i.J 


■ (jS)' 


< e 


(4a) 


-2 ■ • 2 

where e is of order 10 . I" equation (4a), H(«H /u ^ ) is the nondimens ional 

total enthalpy; i is the finite -difference point in the direction along the 

body surface (or s-direction) ; j is the finite-difference point normal to 

the body surface (or n-direction) ; and k is the kth time step. For the 

space-maiching method, convergence criterion of relation (4a) is modified to: 

vf / ..11 1 


maximum 

U 


m __(0 

/3H\* 

W 1-1 


< c 


(4b) 


where Jl is the global iteration count. 

As indicated in the Introduction, the present comparative study of the 
space-^ (due to Davis) auid time-^ (due to MacCormack) marching methods has 
been done by employing these methods as described in references 2 and 3, 
respectively. Since the present methods of solution are identical to those 
of references 2 and 3, these are not presented here. Briefly, since the time- 
marching algorithm^ is explicit and amenable to complete vectorization , it 


has been programed on a vector processor (CYBER 203) . The solutions with 
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the space -marching Inqpllclt-dlffaranoe scheme^ have baan obtalnad on CYBER 
175. Further, in order to spaed up the calculatlcme with the tine-narchlng 
method, the largest possible Courant-Friedrichs-Lei«y (CFL) time-step sise 
for each mesh point has been used In preference to the global minimum CFL 
time-step size. 


DISCUSSION OF RESULTS 

First, a comparison of the predictions made by the space- and time- 
marching methods against the available eiqperimental data is given. Hiere 
are no experimental data available, for the turbulent viscous-shock-layer 
problem (at high Reynolds number) with massive surface blowing considered 
here. For the laminar viscous -shock-layer flow at moderate Reynolds number, 
a comparison between the predictions made by the two methods and the 
eigieriroental data is contained in figure 12 of reference 12. Both the methods 
predict the heating rate for a 45^ half-angle cone which compares favorably 
with the experimental data. For the low Reynolds number case, the predictions 
by the space- and time-marching methods are con^ared against the experimental 
data of reference 13 in figure 2 for a 45° half-angle hyperboloid. Once again, the 
comparison between the predictions by the two methods and the experimental 
data is quite good at Re^»90. At Re^"1035, the predictions depart from the 
experimental data. However, the two predictions compare favorably with 
each other. 

Next, flow-field results and surface heating rates obtained by the 
space- and time-marching methods are presented for the forebody of a 44.25° 
half-angle sphere-cone probe entering the Orton ncmdnal (0.895 + 0.105 He) 
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Jupiter atmosphere under perfect*gas assunption. Results have been obtained 
with and without surface blowing assuming instantaneous transition at the 
first grid point downstream of the stagnation point. The freestream condi- 
tions and other parameters engjloyed in the analysis have been taken from the 

14 

peak heating entry conditions at 51.5 seconds for the 310-kg probe. 

These are given below: 

* 

Nose radius, r„ « 0.352 m 
N 

Free stream Mach number, M =43.76 

OO 

* 

Free stream temperature, T ■ 151.2 K 

00 

* 4 ^ 

Free stream density, « 4,966 x 10 kg/m 

* 

Probe surface temperature, T^ » 4022.80 K 

Ratio of specific heats, y “ 1-4 

* 

Gas constemt, R » 3737.45 J/ka K 

* * * 

Tn p, , 

Reynolds number. Re = — - — ^ « 1.2274 x 10 

Prandtl number, Pr - 0.72 

Mixture molecular weight, m = 2.2246 

Computations by both space- and time-marching methods employed a constant 

* * 

step size As (As /r„) ■ 0.1597 in the direction along the body surface and a 
total of 15 body stations were used. In the direction perpendicular to the 
body, the grid was generated as given by equation (3) and altogether 10 1 point 
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were used with a constant step size An ■ 0.01. Depending on the value of 

6 in equation (3 ) , this step size gives a finer resolution near the body surface 

* * 

in the physical s , n - plane. Two values of 3 (1.1 and 1.01) were 
employed in the computations. For the first value of 3, the first grid 
point off the body surface is located at 0.295 percent of the shoc)c standoff 
distauice as compared with 0.05418 percent for 3 “ 1.01. The ’*arious results 
computed by the two methods are displayed in figures 3 through 9 zmd are 
discussed in the following paragraphs. 

Figures 3(a) and 3(b) contain the surface pressure distributions and the 
shock standoff distances obtained from the two methods with and without 
surface blowing conditions. Except for the region near the tamgency point, 
both these methods give comparable results . Tlie results obtained by the two 
codes at the juncture point should not be taken seriously due to the 
curvature discontinuity there. For an unblown shock layer, the mettiod of 
characteristics does give^^ a fast pressure recovery at the tangency point 
similar to the one obtained by the space-marching method. The fourth-order 
damping employed with the time-marching method smears out the effect of 
surface discontinuity, since a large number of global iterations is used to 
obtain a converged solution. 

Figures 4 through 7 show velocity emd temperature distributions over 
both the spherical (s * 0.64) and the conical flcmk (s » 1.76) portions of 
tlie probe. There is reasonable agreeuient between the various results 
obtained by the two computational methods. However, if the results near 
the wall (figs. 5 and 7) are analyzed carefully, the slopes of the velocity 
and temperature profiles for 3 ■ 1.1 and no surface blowing (m ■ 


0) are 
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seen to be quite different. With the reduction in the value of 0 to 1.01# 
the slopes predicted by the two methods begin to agree quite well. Reducing 
6 further does not result in any noticeable ii^rovement. The surface 
pressure distribution, shock standoff distance, and gross features of the 
velocity and temperature profiles are unaffected by reducing 0 from 1.1 to 
1.01. As will be seen later, the quantities affected most by the mesh 
refinement near the surface are the wall heat transfer and skin friction, 
as expected. Figures 5 and 7 also indicate that with the coarser normal- 
direction mesh (0 >■ 1.1), results obtained by the space -marching method are 
closer to the finer mesh (0 - 1.01) results. With surface blowing, the 
mesh refinement near the surface does not affect the slopes there substan- 
tially. A mesh refinement obtained with 0 > 1.1 appears adequate. Use of 
the largest possible value of 0 without jeopardizing the accuracy of the 
results is important, since the con^utational time is reduced by almost 20 
percent for the time -dependent method when 0 is increased from 1.01 to 1.1 
for the same convergence criterion of equation (4a) . Yet if a large value of 0 
is used with small surface blowing rates, then poor resolution of the 
maximum gradients may occur. 

The eddy viscosity distributions over the spherical and conical 

portions of the probe are given in figs. 8(a) and 8(b) . There is reasonable 

agreement with respect to the maximum values and the general character of 

the distributions obtained by the two methods. It may be mentioned again 

7 8 

that the identical turbulence models ' have been used with the two 


ccxnputational methods . 
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Finally, figs. 9(a) through 9(c) glva the s)iln>frictlonf dimensional 
and nondimensional (Stanton number) )ieatlng rates at the wall. As indicated 
earlier, there is better agreement between the results obtained with a fine 
mesh (6 “ 1.01) for the unblown case. For the blown surface case, the two 
numerical methods give good results for even the coarse mesh (3 ■ 1.1). 
Generally, the results are in better agreement over the conical flanJt portion 
than over the spherical portion or near the juncture point. 

Ccxnputationally , in absence of blowing, the space -marching method is 
quite economical and gives results which compare well with the e]q>erimental 
data. However, with large blowing and arbitrary body shapes, the time-marching 
method has proved to provide solutions without the numerical instability 
experienced with the space -marching method. Without blowing and for analytic 
body shapes, both methods appear equally accurate. 

Further, the time-marching method has been used to predict flows with large 

embedded subsonic regions. Some recent efforts^^ have also been made to 

handle such flows with the space-i.tarching method by making use of certain 

approximations to improve the initial marching profiles. 'Die advantage of the 

space-marching method in such cases lies in tlie fact that for problems 

requiring more computaticmal grid points (either due to the large body size 

or large Reynolds number of the flow) the storage and computational time 

requirements are much less as conpared to the time-marching method. Also, 

with the time -marching method, the ccm^utational time increases substantially 

if more grid points are clustered in the shear layer, whereas with the 

space-marching method, the ccm^utational time is not significantly affected 

by distribution of the grid points in the shear layer. The coa^utaticmal 

time of the time-marching method, however, can be reduced substantially by 

18 


employing its recent implicit analogue. 


11 


Finally, for a viscous raacting flow with nonaquillbrlua chsmlstry, 

19 

there are three tltne->step criteria used with the tlnwmarchlng nethod, 

namely: the Courant-Prledrlchs*'Letfy (CFL) criterion, a viscous dissipation 

criterion, and a species production criterion, therefore, the solution of 

the species equation proceeds with a time step which is different frcm the 

one used with the momentum and energy equations. For the r.cmequi librium 

chemistry calculations, such local differences in time will affect the final 

relaxed solution, since in going through the transient part of the solutions 

one would produce or deplete species at a rate which may not be characteristic 

of the final converged solutions. Unlike the case with equilibrium chemistry, 

the nonequilibrium chemistry cannot correct itself as the velocity and 

temperature fields converge to the steady state values. Ihe space -marching 
20 

method, however, appears to be free from this drawback. 

Ihus, the choice of the numerical method is rather problem dependent. 

Both methods give comparable results for the problems where conqparisons with 
the experimental data have been maae, provided that care is taken in selecting 
the appropriate mesh s^ize near the body suz^'u..e. 

CONCLUDING REMARKS 

Ihe space -marching implicit technique of Davis and the time-marching 
explicit method of MacConnack have been used to analyze the flow field around 
a massively blovm Jupiter entry probe under perfect gas conditions. Hiese 
two techniques have been widely used for the heat shield design calculations 
for such a probe. Certain difference:; *re known to exist in the predicted 
results by tiie two techniques. In order to isolate these differences, a point 
by point con^arison of the various flow-field quantities computed by the two 
schemes has been obtained by using identical computational grids and turbu- 
lence i»)dels. Generally, the results obtained by the two methods are in better 
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•gre«Mnt over the conicel flank portion than over the spherical portion or 
near the juncture point. For the saM decree of accuracy, the space -marching 
scheme takes much less time and would appear to provide accurate results for 
the problems with nonequi librium chemistry, free frcMn the effect of local 
differences In time on the final solution %<hich are inherent in time-marching 
methods. With the time-marching method, however, the solutions are obtained 
for the realistic entry-probe shapes with massive blowing rates, whereas with 
the space-marching technique, it is difficult to obtain converged solutions 
for such flow conditions. Hie choice of the numerical method is, therefore, 
problem dependent. Both the methods give cos^rable results for the problems 
where conqparisons with the experimental data have been m^c, provided that 
care is taken in selecting the appropriate mesh size near the body surface. 
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Nondimensional blowing rate distributions along the surface. 
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Fig. 3 Comparison of time-marching and space-marching solutions both 
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Fig. 4 Comparison of tangentiai velocity profile predictions. 
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Fig. 6 Comparison of temperature profile predictions. 


Tlnve-marching code 


ORIGINAL PAGE 
OF POOR QUALITY 



I 

1 


Fig. 6 Concluded. 
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Fig. 7 Effect of norm^ grid spacing on oredicted temperature profiles. 
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DEVELOP REAL-TIME DOSIMETRY CONCEPTS AND 
INSTRUM£NtAY16>< FOR LONG TE^M MISSIONS 


INTRODUCTION 

Major objectives in the process of developing a rugged portable Instrument 
to evaluate dose and dose equivalent have been achieved. A tissue-equivalent 
proportional counter simulating a 2 micrometer spherical tissue volume has 
operated satisfactorily for over a year. The basic elements of the electronic 
system have been designed and tested. And finally, the most suitable mathe- 
matical technique for evaluating dose equivalent with a portable instrument 
has been selected. Oesigt abrication of a portable prototype, based on 
the previously tested circuits, is underway. 

DETECTOR DEVELOPMENT 

The 5.7 cm detector has been operated for over 15 months with the original 
gas filling. Figure 1 shows the gas gain for a constant anode voltage during 
that time. The long-term drift, less than one percent per month, and short- 
term (les:» chan 24 hour) variations of +^2 percent can easily be controlled by 
adjusting the anode voltage, (see Figure 2). These adjustments will be made 
automatically by a subroutine, illustrated in Figure 3, which compares the 
actual position of a calibration peak with the position of that peak at the 
proper gas gain. The difference is used to calculate a new high-voltage supply 
setting. The actual voltage is provided by a high-voltage supply referenced 
to the output of a digital -to-analog converter. 

The detector will be operated with a combination of gas gain and electronic 
gain which results in a calibration factor of 0.12 keV/ym/channel and a useful 
range of 0.36 to 300 keV/pm. Since the dose mean of the single-event distri- 
bution for X and y rays varies from 0.7 to 3.0 keV/pm, and for neutrons from 
30 to 130 keV/pm, the system is expected to detect photon as well as neutron 
dose. However, the stainless steel vacuum chamber surrounding the detector 
will distort the response as a function of photon energy. Also, the photon- 
induced events cannot be distinguished from very low-energy neutron events. 
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FIGURE 1 . Detector Gas Gain at A Constant Anode Voltage as 
a Function of Time After Sealing the Detector 

Thus, though the system is expected to accurately measure dose due to photons 
and low-energy neutrons, estimates of the mean quality factor become less 
certain when low-energy neutrons are abundant. 

Specifications for the second detector, intended to detect high-energy 
heavy particles, have been completed. It will use the same basic design and 
materials but will be 12.7 cm in diameter in order to provide approximately 
five times the counting rate. 

ELECTRONIC SYSTEM 

The electronics for a system using two detectors is outlined in Figure 4. 
Each detector is supported by a dedicated high-voltage supply, amplifiers and 
analog-to-digital converters (ADC's), but the two detectors share a multichannel 
analyzer (MCA) and microcomputer. The ADC's output is used as a MCA memory 
address. Each event results in a "one" being added to the memory content at 
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FIGURE 2 . Gas Gain Versus Anode Voltage 


the corresponding address. As illustrated in Figure 5, the entire MCA content 
can be periodically transferred to the microcomputer memory. The dosimetric 
quantities can then be computed without interrupting further data collection. 
An absolute time clock in the microcomputer is used to initiate the dose cal- 
culation, data storage, calibration and other functions. A CMOS successive 
approximation ADC chip and other CMOS components are used in the ADC circuit, 
in Figure 6, to minimize power consumption. 
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